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FOREWORD 

This  is  Volume  II  of  the  final  report  under  Contract  No.  Nonr  5006(00). 
This  report  is  issued  in  order  to  present  the  linear  theory  (in  its  most 
general  form)  of  water  waves  produced  by  explosions  and  to  specify  the 
assumptions  which  prevail  in  performing  the  Mono  Lake  deep  water 
predictions.  The  most  general  form  of  the  linear  theory  has  never 
..-appeared  in  the  literature  prior  to  this  report  and  it  is  evident  that  it 
has  not  been  employed  to  its  fullest  extent. 

This  report  was  prepared  under  the  direction  of  Dr.  B.  Le  Mdhaui£, 
Associate  Director  and  Project  Manager  for  this  program.  Miss  Lois 
Webb  contributed  to  the  development  of  the  computing  techniques  and 
the  programming  was  accomplished  by  Mr.  H.  Lee  Butler. 


ABSTRACT 


The  objective  of  this  project  was  to  predict  the  water  waves  generated 
by  the  Mono  Lake  field  tests.  These  tests  consisted  of  detonating  ten 
TNT  explosions  of  approximately  9.  500  pounds  each,  and  measuring 
the  water  waves  and  run-up  produced  by  these  detonations. 

Volume  X  of  the  final  report  contains  these  predictions  and  a  brief 
description  of  the  theoretical  and  empirical  methods  employed  in  per¬ 
forming  the  predictions.  This  volume  presents  the  general  form  of  the 
linear  theory  of  water  waves  produced  by  explosions  in  water  of  con¬ 
stant  depth.  The  theory  is  developed  in  detail  so  that  the  full  extent 
of  the  assumptions  made  in  performing  the  predictions  for  Mono  Lake 
can  be  realized.  It  is  also  evident  that  several  natural  extensions  to 
the  presently  employed  theory  can  be  made  which  will  result  in  more 
realistic  initial  conditions  and  which  should  present  a  more  reliable 
theory.  Examples  of  the  various  solutions  that  depict  the  differences 
between  a  symmetric  surface  deformation,  a  symmetric  time-dependent 
surface  deformation,  and  an  asymmetric  time-dependent  surface 
deformation  are  given. 
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i.  INTRODUCTION 


The  linear  theory  of  water  waves  produced  by  explosions  has-been 
utilized  by  many  investigators  (Lamb  1945;  Van  Dorn  1959*  1963; 
Kranzer  and  Keller  1955;  Prins  1956;  Whalin  1965;  Kajiura  1963; 

Unoki  and  Nakano  1953;  Stoker  1947;  Penney  1950;  Fuchs  1952; 
Kirkwood  and  Seeger  1950)  during  the  last  20  years.  The  basic 
mathematical  foundation  for  the  generalized  linear  theory  has  been  in 
existence  for  a  longer  period  of  time.  It  involves  the  use  of  a  time- 
dependent  Green's  Function  which  satisfies  the  Laplace  equation  with 
the  appropriate  linearized  boundary  conditions  and  the  recognised 
assumptions  that  linearization  involves.  However,  this  theory  in  its 
general  form  has  escaped  application  with  one  notable  exception 
(Kajiura.  1963).  who  applied  the  theory  to  calculate  the  leading  wave 
of  a  tsunami  and  gave  various  examples  of  the  theory  to  show  that 
certain  special  cases  reduced  to  the  well  known  solutions  (Kranzer 
and  Keller.  1955)  for  explosion  waves. 

This  report  formulates  the  general  theory  and  applies  it  to  water  waves 
produced  by  explosions.  The  asymptotic  solutions  are  developed, 
along  with  examples,  for  several  special  cases  and  the  full  impact  of 
the  use  of  the  time-dependent  Green's  Function  is  thoroughly  discussed. 
The  advantage  of  using  this  theory  in  the  case  of  one  detonation  and 
multiple  detonations  is  pointed  out. 

The  axially  symmetric  water  craters  presently  being  used  to  predict 
explosion  waves  many  times  bear  little  resemblance  to  the  actual 
craters  created  (Van  Dorn  1963,  Whalin  1965);  these  craters  are 
generally  termed  "pseudo"  craters.  The  use  of  a  time-dependent 
initial  deformation  will  undoubtedly  result  in  much  more  realistic 
crater  parameters.  These  parameters  should  be  related  to  weapon 
yield  and  location  as  has  been  accomplished  for  stationary  symmetric 
water  craters  (Whalin,  1965).  The  theory  in  its  most  general  form  can 
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be  applied  to  predict  the  wave  spectrum  generated  by  an  exotic  shaped 
source  resulting  from  multiple  detonations  which  may  be  employed  to 
produce  wave  amplification  at  a  certain  point  or  in  a  specified  area. 

The  conclusion  of  this  report  is  that  the  prediction  methods  presently 
employed  and  used  for  the  Mono  Lake  program  should  be  modified  to 
account  for  a  time- dependent  source  of  disturbance.  Further,  the 
empirical  relationships  between  the  cavity  parameters  and  the  weapon 
size  and  location  should  be  established  as  they  have  been  for  a 
stationary  initial  deformation.  Furthermore,  one  can  probably 
establish  a  model,  which  will  account  for  the  evaporation  of  water  as 
well  as  the  subsequent  return  of  water  from  the  plume  collapse.  In 
addition,  the  general  solution  for  an  arbitrary  shaped  initial  deforma¬ 
tion  can  be  used  to  predict  the  wave  spectrum  generated  by  multiple 
explosions. 


2.  THEORETICAL,  FORMULATION 

The  mathematical  medal  employed  in  the  generalised  linear  theory  of 
water  waves  generated  by  explosions  is  based  on  the  linear  theory  of 
surface  waves  in  an  inviscid  incompressible  fluid  of  constant  depth. 

This  model  allows  the  source  to  be  asymmetrical  and  time -dependent 
in  its  most  general  form.  Consequently,  a  mathematical  model  is 
available  in  which  the  actual  wave  generation  mechanism  can  be 
simulated  much  more  precisely  than  that  generally  used  in  present 
applications. 

The  water  is  assumed  incompressible,  of  constant  depth  d;  the  origin 
of  the  Cartesian  coordinate  system  (x'»  y',  z ' )  is  at  the  undisturbed 
free  surface,  the  vertical  axis  z*  is  positive  upwards,  and  the 
variables  are  transformed  to  dimensionless  form  (i.  e. ,  x  =  x//d, 
y  =  y'/d,  t  =  t'  Vg/d,  V  =  tj'/d,  V  =  vVv  gd  ;  <0  =  o'/d  VgdT 
P  =  (pVpl/gd;  etc.). 

The  kinematic  and  dynamic  conditions  at  the  free  surface  are 

<PZ  *  *  =  0  (1) 

<0t  =  -  tj  -  p,  a  S  0  (2) 

and  the  bottom  condition  is 

pz  a  WB.  8  a  - 1  (3) 

where  Wg  is  the  assumed  bottom  velocity  corresponding  to  a  bottom 
deformation.  The  above  conditions  are  subject  to  the  restrictions  of 
assuming  a  linear  approximation.  That  is,  the  deformation  at  the 
surface  and  bottom  are  assumed  small  compared  with  the  wave  length, 
A/,  and  water  depth,  d,  and  the  condition  s/X/6/dJ  <  1*. 
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The  time-dependent  Green's  Function,  G,  is  a  solution  of  the  Laplace 
equation 

V2G  =  0,  0  >  z  >  -1  .  t  *  r  (4) 

satisfying  the  free  surface  condition 

G_  +  G  =  0,  z  =  0  (5) 

tt  2 


and  the  bottom  condition 

G  =  0,  z  =  -1  (6) 

s 

It  is  required  that  G,  G  ,  G  ,  G  ,  G. ,  G.  ,  G.  »  and  G.  be  uniformly 

x  y  z  t  tx  ty  tz 

bounded  for  every  t  at  x, y  *• m.  In  addition,  (G  -  1/R)  must  be  bounded 
at  (xo,  yQ,  zQ)  where 

R2  =  (x  -  xo)2  +  (y  -  yo)2  +  <z  -  zQ)2  (7J 

At  time  t  -  T  the  following  initial  conditions  are  assumed 


G  =  Gt  =  0,  *  =  0 


The  above  equation  and  conditions  are  sufficient  to  determine  G  uniquely. 
From  Stoker,  the  Green's  Function  for  the  case  of  three-dimensional 
motion  in  water  of  finite  depth  is 


C(W  z0;Tlx’  Y*a» 


p®  j 

*>  =  I  FoTh'm  [»i**m{l  -  |z-zol}-sinhm{l  +  (z  +  zo)} 


+  {l  -  coshy(t-T)}-££^~-j£  cosh  m{l  +  z)  coshm(l  +  zo)]  dm, 

0  >  z,  z>  - 1 
o 
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where 


(*  -  *e)2  +  (y  -  y0)2 


00) 


and 


y  -  m  tanh  tn 


(ID 


Note-that  G  is  symmetric  with  respect  to  (x  ,  y  ,  z  )  and  (x,  y,  s)  and 

o  o  o 

also  t  and  r.  Therefore 


o' '  o 


(12) 


From  Green's  formula  and  the  above  Green's  function 
<PrU.y.«;r)  =  i  Jf  (op  -  <PrG  )  aso 

o  O  O  * 

O 

•  Lf(c®r.  '  fT°z  )  ,dso  <13> 

5  o  o  Z  *-l 

o 

where  dSQ  =  dxQ  dyQ  and  the  integral  on  the  lateral  boundary  in  water 
vanishes  because  of  the  condition  imposed  on  G. 

The  integration  of  Eq.  13  with  respect  to  T  from  0  to  t  and  substitution 
of  conditions  in  Eqs.  1,  2,  5,  6  and  8  for  ip  and  G  yield 

«<*.  y.  «.t) .0)  •  .  - Orn)  +  f  KV#T]  «SD 

V 

■h!JS0\°\r)E  c.  ,*••>=<.  ««> 

o 

5 


r 


The  first  integral  in  Eq.  14  represents  the  contribution  c£  th®  curias® 
condition  and  the  latter  integral,  that  o£  the  bottom  cessation. 

At  the  surface  z  =  0,  Eq.  9  becomes 

i  (*•  2  »  J  (m3?) 

G(xc,yo,  z0;r | x, y,  0, t)  =  -j  {l  -  cos y(t-r)}m cosh m(l  +  so)^--m  dm 


(15) 


and  since  the  surface  elevation,  1},  is  given  by  Eq.  2,  the  substitution 
of  Eq.  14  into  Eq.  2  yields 


”  W  4K+I2  +  I3|d3 


(16) 


where 


U  * 


I,  ■ 


I*  * 


(Gt<VCrt’1U:  *°  *  °'  "  °’ 

J0  P  °m  iT  *  " 

J*  pr  crt iT  '  (*,OTt)T=0!  *o  *  *•  *  "  *•  •”d 

j' G.  \r  ir-  01 

-£0‘rVr'(°‘ %),.<>’  ‘o'-'-"0- 


(1?) 

(18) 

(19) 

(20) 
(21) 


Note  that  Eq.  17  represents  the  contribution  of  the  initial  velocity  and 
elevation  of  the  water  at  the  surface,  Eq.  18  or  Eq.  19  the  contribution 
of  the  surface  pressure,  and  Eq.  20  or  Eq.  21  the  contribution  from  the 
bottom  deformation. 
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3.  ASYMPTOTIC  SOLUTIONS 


This  section  presents  the  asymptotic  solutions  obtained  by  an'  spplicstica 
of  the  method  of  stationary  phase  to  evaluate  the  integral  in  Eq.  16. 

The  special  cases  for  which  the  asymptotic  solution  is  developed  are 
a  symmetric  initial  impulse*  a  symmetric  initial  deformation*  a 
symmetric  time-dependent  initial  deformation  and  an  asymmetric  time- 
dependent  initial  deformation.  In  each  case  the  source  is  given  as  an 
-arbitrary  function  where  special  cases  are  easily  investigated.  Examples 
of  a  special  case  appear  for  each  assumed  initial  condition.  The 
condition  of  a  symmetric  initial  impulse  (Section  3. 1)  and  a  symmetric 
initial  deformation  (Section  3. 2)  were  solved  by  Kranzer  and  Keller 
(1955)  and  are  reproduced  here  for  completeness.  However,  the 
developments  of  Sections  3.3  and  3.4  are  original  and  bear  further 
investigation  as  far  as  application  to  actually  performing  predictions 
are  concerned. 


3. 1  Axially  Symmetric  Initial  Impulse 


The  special  case  of  an  axially  symmetric  initial  impulse  was  solved  by 
Kranzer  and  Keller  by  the  method  of  integral  transforms.  When  the 
proper  assumptions  are  employed  their  solution  is  obtained  from  Eq.  16. 
Cylindrical  coordinates  are  used  and  the  initial  conditions  become: 


a)  The  water  surface  is  initially  undisturbed 


%  (ro*  9o’ T)  *  0 


The  initial  velocity  at  the  surface  at  time  t>  0  is. 

s  «  0 
o 

»  *  0  *  0 
°  T  *  t  «  0 


Consequently  I  j  «  (c^  - 


T=0;j  *0,B*0 

o 


7 


The  applied  pressure  at  the  surface  is  isipsLsiv®  a»  7  *  C  , 

i.  e. ,  p  =  I  £{7);  the  refers 

o 

*2  *  Jjp.Gm iT  *  0>  gtt)  T=t-  *„  ■  °«  *  *  0 


■  W  °rrf  r  •  0.  .  .  0.  •  0 


-  d)  1 3  *  0 

Upon  substitution  in  Eq.  16  and  expressing  the  wave  amplitude  in 
cylindrical  coordinates 

®  *2lf  *  • 

ti(t,  0,  t)  =  ^  J  2y  m  sin  yt  |  J  J  IQ(ro)  JQ(mr)  rodrod0o}-dm  (22) 


or  upon  integration  with  respect  to  tq  and  0q  the  wave  amplitude  becomes 

«w  _ * 

Tj(r,  0,  t)  =  -  J  m  Vm  tanh  m'  !T  (m)JQ(mr)  sin  (Vm  tanh  m  tjdm  (23) 


where 


T"(m)  =  f  I  (r  )  J  (mr  )  r  dr 

O  «/•  OOO  OOG 


The  above  integral  is  evaluated  by  the  method  of  stationary  phase  upon 
the  assumption  that  the  frequency  of  oscillation  of  J o(mr)  and  the 
circular  function  is  large  compared  to  that  of  I^(m).  Then  Jo(mr)  is 
approximated  by  an  asymptotic  expansion  and  upon  integration 

»j(r,t)  «  -2— -  "\J  sin  (mr  •  Vm  tanh  tn  t)  (25) 


where 


^m) 


nh  m 


2  cosh  m 


tanh  rn 


B 


Table  I  in  Section  3. 2  gives  T^{m)  for  various  initial  impulses  I0(?0). 
Figure  1  shews  an  example  c£  the  wave  amplitude  as  a  function  cl 
time  for  the  initial  impulse  in  example  2  of  Table  I  at  cue  observation 

station. 

3.2  Axially  Symmetric  Surface  Deformation 

The  asymptotic  solution  for  an  axially  symmetric  surface  deformation 
was  solved  by  Kranser  and  Keller  for  several  special  cases.  However* 
the  theory  is  presented  in  the  following  pages  from  the  more  general 
mathematical  formulation  given  in  Section  2.  When  cylindrical 
coordinates  are  used*  the  initial  conditions  become: 


a) 

b) 


c) 


The  water  surface  is  initially  undisturbed  and  is  defined 

by  the  function  fj  (r  ). 

o  o 

The  initial  velocity  of  the  surface  at  time  t  =  0  is 


s  s  0 
o 

z  =  0  =0 


1 r  *  t  =  o 
*1  ”  *  Grt  Vro^T=0,z  =0,  z=0 


=  [-  J  2  cos  (yt)  m  J0(mr)  dm]  V0(rQ) 


d) 


*  1  3  *  0 


Upon  substitution  in  Eq.  16  and  expressing  the  wave  amplitude  in 
cylindrical  coordinates 

Vir.t)  -  ^  Jq  2m  cooyt{£  £  n^tymr^dr^ejdm 


(27) 


or 


=  i?J0  2mcoa  ffJ0  ’}o(ro){Jo(mr)Jo(tnro)  +  1Jn(m?)Jn(ciro,c03  ^°'do}} 

r  dr  d9  1  dm 

O  O  OJ 

(28) 

however 

J 

|  cos  n  (6  -  0  )  d0A  =  0 

0  o  o 

(29) 

therefore 

* r.  t)  =  -jq  m  Jo(mr)  cos  yt  [£  » ?0(ro>  J0(mrQ)  rQ  drj  dm 

(30) 

where 

J 

F0  V'o1  ,0<mro)  'o  drO  *  %'”> 

(31) 

which  ia  the  aero  order  Hankel  Transform  of  n  (r  ).  The  wave 

o  o 

amplitude  becomea 


lj(r,  t)  »  -  f 


rjo(m)  m  cos  (j\f 


m  tanh  m  t)  JQ(mr)  dm 


(32) 


The  above  integral  ia  evaluated  by  the  method  of  stationary  phaae  where 
it  ia  assumed  that  »J0Cm)  has  a  much  lower  frequency  of  .oscillation 
than  JQ(mr)  and  the  circular  function.  This  assumption  is  a  valid  one 
for  the  types  of  initial  deformations  which  are  representative  of  under* 
water  explosions,  as  can  be  observed  from  the  examples  in  Table  X. 
JQ(mr)  is  approximated  by  an  asymptotic  expansion  for  large  argument 
and  the  integral  becomes 


»K*.t) 


m 


T?o(m)  cos  (mr  -  tanh  m  t) 


(33) 


U 


where 


#m)  ■ 


1  -ftsssh  m  .  1  /'  ra 

1  V  m  Zcoih^m 


r 

r 


(34) 


The  following  table  gives  T]  (m)  for  various  initial  deformations  t)  {r  ). 

O  O  O 

The  wave  amplitude  may  be  computed  either  as  a  function  of  time  with 
rconstant  or  as  a  function  of  r  with  time  constant. 

#/ 

Figure  2  shows  the  wave  amplitude  as  a  function  of  time  for  example  3 
of  Table  I  at  one  observation  station. 


3.3  Axially  Symmetric  Time-De pendant  Surface  Deformation 


There  has  been  little  investigation  into  the  effect  of  using  a  time- 
dependent  symmetric  surface  deformation.  This  section  presents  the 
theory  with  several  specific  examples  given  in  detail.  The  examples 
can  be  compared  with  those  where  the  source  is  independent  of  time. 
The  intial  conditions  become: 


a) 


The  water  surface  is  initially  undisturbed  and  is  defined 
by 


VvT)*{ 


Vro*  ,in  I  7?  * 

0, 


0  *  f  *  T* 
T  >  T* 


b) 


It  is  assumed  that  the  variables  r  and  r  are  separable  and 

o 

r*  should  be  taken  as  the  dimensionless  period  of  the  first 
expansion  of  the  bubble  or  the  water  cavity  for  the  case  of 
a  surface  explosion. 


The  Initial  velocity  of  the  surface  at  time  r*01s 


•  0 
-  0 

»  t «  p 


dWT) 

a*r 


T»0 


iff  W 
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TABLE  I 

HANKEL  TRANSFORM  OF  VARIOUS  INITIAL  DEFORMATIONS 


t?Jr  )  I  (rQ) 
o  o  _ o  o 

~fj -  or-T — 

o  o 


%(m)  "Um> 

-rT~  or  ”1 

'o  o 


1.  1,  r  s  R 

o 

0,  *„>  R 


[‘  •  t  (tt)2], 


.  ro  >  V2R 


4  „  2 

-r_  4r 


+  1^r°sVlR 


.  r  >  VIr 

o 


r  r  2i  *3/2 
4.  [l  +  2  ] 


£  I,MI 


— r  J.  (  t/ZtoR) 
m2  2 


J4(  V-3  mR) 


R2  .-RmVT 

T  c 


5.  e 


V~2  r/R 


2(1  +  R2m2/2)” 


6.  e-<r/R> 


(?) 


.  VgJ  .. 


c) 


'l  *  (°t«v) 


and 


T=0,  2  =0,2=0 

o 


G.  =  j  \  sin  yt  mJ  (mr)  dm 
*  -0  7  ° 

d)  x2.i3.o 

Upon  substitution  in  Eq.  16  and  expressing  the  wave  amplitude  in 
cylindrical  coordinates 

=  ■2ayB**EI*  J0Jote',Vro,(l^)rodro d9o] *“ 

*  7F*  J0  m*y  **  [£  J0  {,<,<mrlJ„(mroH2iif{n(mrlJ»<,"ro,co» 

ri(r)rdrd0ldm 
o  o  o  o  oJ 


=  +  Tr*  J0  -^r—  Jo(mr)  [  J0  Vro>  Jo(mro)  ro  dro]dm  (35> 

However,  the  integral  with  respect  to  r  is  merely  the  zero  order  Hankel 

c 

Transform  of  the  initial  deformation  TJo(rQ)  will  be  denoted  by  tj^(m). 

.  -®  mJ  (mr)  _  _ 

Vir,  t)  =  ——z  -  ■  - ; . .  n  (m)  sin  ( Vm  tanh  m  t)  dm  (36) 

ZT  0  -^m  tanh  n  ° 

Upon  evaluation  by  the  method  of  stationary  phase 

*Kr»t)  =  — ■°—-  ,in  (mT  "  V m  tanh  m  t)  (37) 

r  V  tanh  m  “  ' 

where  0  and  <p'  are  defined  as  before.  Examples  of  rj£(m)  are  found  in 
Table  I. 
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Figure  3  shews  the  wave  amplitude  as  a  function  of  time  at  csa 
observation  station  where  is  given  by  example  3  cf  Tcbls  I. 

3.4  Asymmetric  Time-Dsner-dant  Surface  Deformation 

In  addition  to  assuming  the  s  ource  of  disturbance  to  be  time-dependent, 
it  is  relatively  straightforward  to  avoid  the  assumption  of  axial 
symmetry.  However,  in  the  case  of  explosion  waves,  this  assumption 
is  not  necessary  unless  an  unsymmetric  cavity  is  produced  by  an  in¬ 
complete  detonation  or  by  multiple  detonations.  The  initial  condition 
becomes: 


The  water  surface  is  initially  undisturbed  and  is  defined  by 

A(r^8(0o,T(T)* 

1o(ro'eolt)  =  L 


0  £  T  <  T*, 

0  S  8  <  2» 

t  >r,° 

*o>R 

b) 


It  is  assumed  that  the  variables  are  separable  and  that  r* 
is  a  parameter  which  should  be  related  to  the  time  of  the 
initial  expansion  of  the  bubble. 

The  initial  velocity  of  the  surface  at  time  T  =  0  is 


<f> 


z  =  0  d  T)  (r  ,  8  ,  T) 
o  _  o  o  o 


z  =  0 
T  =  0 


"HF 


T=0 


c) 


1  1  *  GtV 


*  Vro>  Tr(0) 


and 


T=0,  z  *0,  z=0 
o 


08 

G  *  J  |  sin  y  t  m  J  (air)  dm 

•  7  O 


d) 


W° 
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Upon  substitution  in  £q.  16  and  expressing  the  wave  amplitude  in 
cylindrical  coordinates 

«*•••«  [^{j.tarl  *.<«»,,)♦  2|  JB(rr,t)Jn(n,ro)cc3  riS-S,,)} 

•  T?o(ro)  e{0o)rodrod9Q]dm  (33) 

Upon  performing  the  integrations  with  respect  to  and  0^ 


T  (0)  — 

rtr,  0,  t)  =  J^i!l^l21[jo(mr)^(m)  I#  +  2  ^Jn(mr)l^(m)  ln(0)]dm  (39) 

re 

v  (m)  =  I  f)  (r  )  J  (mr  )  r  dr  *  aero  order  Hankel  transform 
J0  °  °  °  000  of  n  (r  > 


where 


'o  o' 


Vm)s  J0Vro,Jn(mro,rodro  “  order  Hankel  transform  of 


Vro> 


.  f 


Jq  6<V  d9o 

f  A  ffsinn0 

I  (©)  *  {  n 

n  VB  #  cos  n0 

n 


(40) 

(41) 


and  the  A  are  coefficients  of  the  Fourier  sine  series 
n  • 

representing  0(0o),  i.  e. ,  if  0{0Q)  ^p  **n  P®0  or  ** 

0(0  )  s  B  cos  p8  .  then  the  B_  are  the  coefficients  of 
o  p=o  p  °  n 

the  Fourier  cosine  series  and  p>n, 

The  above  integral  is  evaluated  by  the  method  of  stationary  phase  and 

»J(r,  0,  t)  «  -  TT(m)  +  2  f  tT  (m)  I  (0)} 

2ir  r  Vtanh  *"  *  -P  (m)  100  ,  n  n  J 

f  n*  1 

sin  (mr  •  ■ym  tank  m  t) 

Ihe  following  table  illustrates  the  various  parameters  in  the  above 
solution  for  specific  initial  deformations. 


(42) 
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Figures  4,  5,  6,  7,  8,  9,  10,  11,  end  12  chew  file  wave  amplitude 

as  a  function  cf  time  at  sis  observation  stations  which  are  equidistant, 

r,  from  the  origin  and  at  various  locations  0  ca  the  surface  where  the 

initial  conditions  are  taken  from  example  3  cf  Table  H.  A  comparison 

with  Fig.  3  of  the  previous  section  shows  the  effect  of  considering  an 

asymmetric  time  dependent  source  where  the  only  difference  is  the 

function  8(8  ).  It  should  be  noted  that  the  consideration  cf  an  aeym* 
o 

etric  source  is  of  little  value  when  considering  the  effects  of  a  single 
—detonation  unless  it  happens  to  be  an  incomplete  and  asymmetrical 
detonation.  However,  the  method  is  certainly  applicable  to  multiple 
detonations. 

The  physical  significance  of  considering  the  initial  deformation  cf  this 
form  should  be  clearly  understood.  It  appears,  at  a  cursory  glance, 
that  we  have  obtained  an  asymptotic  solution,  valid  at  large  distances 
from  the  source,  which  yields  the  solution  for  a  deformation  that  forms 
in  a  prescribed  manner  from  time  0  to  T*,  and  that  the  solution  is 
valid  for  all  t  a  r*.  However,  this  is  not  precisely  the  case.  Actually, 
we  are  solving  an  initial  value  problem  where  we  prescribe  the  shape 
of  the  water  surface  and  the  velocity  field  at  the  surface  at  time  t  *  T  *  0. 
The  velocity  field  is  obtained  from  the  prescribed  time  dependent 
deformation  from  time  0  to  T*.  The  result  is  that  the  solution  obtained 
will  be  a  solution  for  a  class  of  initial  conditions,  precisely  that  class 
which  has  the  same  surface  e.iape  and  velocity  field  at  the  surface  at 
time  T  *  0.  This  class  is  described  by  the  set  S  consisting  of  elements 
T(r)  such  that 

T(0)  «  0  and  Tt<0>  - 
Examples  of  this  set  include 

T(t)  ■  sin  ^  +  Tn  +  rn  cos  ar,  n  *  2 

-sin|  £*TnsbTP,  n*2 

_  w  r 

“  I  T* 
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That  is,  the  solution  does  not  depend  on  the  precise  manner  in  trlsich 
the  deformation  builds  up  from  time  0  to  T*  but  only  ca  kcw  it  starts 
to  build  up.  Furthermore*  if  the  solution  were  evaluated  at  t  *  t*»  it 
would  not  yield  the  precise  deformation  assumed  (or  thought  to  be 
assumed).  However*  it  can  be  expected  that  the  solution  probably 
would  not  deviate  significantly  from  T}Q(ro,  since  the  assumed 

initial  conditions  are  physically  reasonable  conditions  and  the 
mathematical  model  also  is  physically  reasonable. 


4.  METHOD  OF  SOLUTION  IN  THE  AREA  NEAR  THE  SOURCE 


The  asymptotic  solutions  (method  of  stationary  phase)  presented  in  the 
previous  section  are  not  valid  in  the  area  near  the  explosion.  The 
integral  with  respect  to  m  must  be  evaluated  in  some  other  approxi¬ 
mate  manner  and  the  method  proposed  by  Whalin(1965)  for  a 
symmetric  water  cavity  is  also  applicable  with  only  slight  modification. 
The  equation  to  be  solved,  in  general,  is 

tj(r,  8,t)  *  ~  JQm  y°2:t[1oJo(mr)  %{m)  +  Z\  Jn(8)Jn(mr)^(m)]dm 

where  the  quantities  I  ,  “(m),  1  (8),  and  y  are  the  same  as 

o  o  n  n 

previously  defined.  The  proposed  method  of  evaluation  is  to  fit  a 
fourth  degree  polynomical  to  each  half  oscillation  of  the  integrand  and 
integrate  directly.  This  method  requires  that  all  computations  be 
performed  in  double  precision  and  the  upper  limit  of  integration  may  be 
chosen  to  bound  the  error  to  almost  any  predetermined  limit. 


The  above  integral  is  to  be  approximated  in  the  following  manner 

tj(r,  8,  t)  »  ^  J  g(m)  dm  +  J  “"'J'Jz  +  ^  J  k(m)  8m 


where 

a  ■  a  small  number  such  that  the  oscillatory  functions  of  the 
integrand  may  be  approximated  by  asymptotic  expansions 
near  m  •  0. 

8  *  a  value  of  m  large  enough  so  that  the  contribution  of  the 
latter  integral  is  negligible  and  in  addition  is  a  aero  of  the 
integrand,  This  must  be  carefully  chosen  and  depends  upon 
tj£(m)  and  tjT(r.j). 

g(m)  ■  the  Integrand  when  approximated  by  asymptotic  expansions 
for  small  m 
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Having  chosen  the  proper  limits  a  and  0,  the  first  integral  may  be 

integrated  directly  and  the  latter  integral  is  negligible.  The  second 

integral  will  be  approximated  by  a  polynomical  fit.  nfj(m),  over  each 

half  oscillation  of  the  integrand. 

3/2  f  £  % 

nf(m)  «  - - 5£*LH{ioJo(mr)  TJo(m)  +  2^  IJ6)  Jn(mr)  r^(m)} 

where 


of(m)  ■  Jo(mr)  T?o(m),  n  «  0 


nf(m)  ■  2  In(8)  Jn(mr)  tjjm),  n  *  1 
n? 


,flm)  -  {.A' m4  ♦  „B/m3  ♦  4  . 


a  <  m  <  i  *  1 
n  i 

m.  ,  s  m  <  m.,  i  ■  2,  3, ...  1 
n  ni  n 


where  the  m  and  the  ordered  seros  of  the  integrand,  f(m),  and  I  is 
n  *  n  n 

the  number  of  zeros  of  the  integrand  from  a  up  to  and  including  0. 

Upon  substitution  and  subsequent  integration,  the  second  integral  and 


consequently  the  wave  amplitude  becomes 

I P  CaJ/  a/.  Q/}\  -®t/  7/ 


*.•/.  VW'2  -  *,,2)t5r(nm!/2-  •  *!/2) 

«•  nU:W/2-«3/2)*^(X/2-*1/2) 

•  —  1  "I  ^  ^ 

p(  gh)7“  .  WnAl/  9/2  9/2\.nBi/  7/2  7/2N  +  nCi/  3/2  3/2 

i^LT'\nmi  nml- 1  /  +  "7"  \nml  *  nmi- J +“T  Vnmi  *nmi- 1 


!♦■£(  roJ«.  „?«)♦  E,'(  m‘/2-  ">!'?)] 

3  \n  1  n  1-1/  nl\ni  n  i-l/J 
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where 


1 

I 


nAl  *  oAl 


oBi  "  *  (nml  +  nmi-l)  nAi  +  nBl 

nC'i  '  nmi  nmi-l  aAi  *  (nmi  +  nral-l)  nBi  +  nCl 

vPl  *  nmi  nml- 1  nBi  *  (nmi  +  n“i-l)  nCl 


nE4  “  nmi  nmi- 1  nEi 


,  -64  n  ,/nn’i  +  3nmi-n 

V~.  _  i-,m  n‘{ — i - r.H- — i — i 


m,  +  m. 


n  *  (  m.  -  m.  .) 

n  i  n  i- 1 


-  ,3  m.  +  m. 

♦!.«(  n--‘-T-°--)] 


»Bi 


64 _  ^5nml  **  f  ^nml  T  |  nml- 1  ^ 


m.  +  3  _m. 


<iml  -  nnVl> 


m.  +  m. 


/  , /n“‘i  n‘i-n 

\nmi  +  nmi- 1  n  \ - 2 - / 


,3  m.  +  5  m. 


.3  m,  +  m. 


nCi 
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Example*  of  this  integration  method  are  shown  in  Fig*.  13  to  IS  for  the 
case  of  a  parabolic  impulse  and  further  details  can  be  found  in  Whalin 
(1965).  The  method  is  relatively  easy  to  use  and  the  accuracy  of  the 
integration  procedure  is  well  established. 
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A  Comparison  of  the  Asymptotic  Solution  for  a  Parabolic 
Impulse  With  an  Alternate  Integration  Method 


5.  CONCLUSIONS 


The  preceeding  theoretical  developments  present  the  limit  to  which  one 
may  employ  a  linear  constant  depth  model  to  the  problem  of  water 
waves  generated  by  explosions.  The  present  state  of  the  art  is  at  a 
point  where  much  can  be  gained  by  the  natural  extensions  of  the  theory 
presented  in  Sections  2,  3,  and  4. 

«*- 

Data  and  methods  are  available  which  will  enable  one  to  relate  the 
parameters  of  a  time-dependent  symmetrical  cavity  to  the  sise  and 
location  of  the  weapon.  In  addition,  there  are  sufficient  data  available 
from  the  Mono  Lake  tests,  Hydra,  and  other  WES  tests  to  analyse  the 
usefulness  of  considering  cavity  parameters  which  result  from  this 
model.  Figures  2  and  3  indicate  the  difference  obtained  in  the  wave 
train  by  considering  the  time  dependence.  This  order  of  magnitude 
difference  in  the  envelope  amplitude  is  a  result  of  the  fact  that  a 
symmetric  stationary  cavity  contains  only  potential  energy  and  the 
assumed  time-dependent  cavity  from  time  0  to  r*  will  contain  approx¬ 
imately  the  same  potential  energy  at  Mme  T*.  However,  in  addition, 
the  fluid  contains  kinetic  energy  at  this  time  and  henc*  the  resulting 
vast  difference  in  the  envelope  amplitudes.  A  thorough  investigation 
of  this  form  of  initial  conditions  should  prove  beneficial.  Also,  it  has 
been  pointed  out  by  Van  Dorn  that  one  could  probably  even  consider  the 
effect  of  the  plume  collapse  by  the  proper  choice  of  initial  conditions 
and  skillful  use  of  land  cratering  data. 

The  utility  of  considering  an  asymmetric  time-dependent  form  as  an 
initial  condition  is  obvious.  If  the  effects  of  multiple  explosions  are  to 
be  studied,  then  this  is  the  only  valid  method.  One  can  not,  in  general, 
merely  add  the  solutions  of  several  symmetric  surface  deformations 
after  obtaining  the  asymptotic  solution  for  each.  The  total  deformation 
of  the  source  area  must  be  considered  initially,  then  transformed  and 
Integrated  approximately  by  the  method  of  stationary  phase  (m.  e.p)  to 
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obtain  the  correct  asymptotic  solution.  It  also  should  be  noted  that 
tiie  asymmetric  model  is  applicable  to  the  predictions  cf  the  wave 
spectrum  generated  by  an  earthquake  when  applied  to  circumstances 
consistent  with  the  assumptions  of  the  mathematical  model. 

Predictions  at  Mono  Lake  were  made  from  example  3  of  Table  I  which 
is  a  symmetric  surface  deformation  surrounded  by  a  lip  such  that  the 
volume  contained  in  the  lip  is  equivalent  to  the  volume  of  the  cavity, 
i.  e. .  no  net  addition  or  removal  of  fluid  from  the  system. 
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NO  1  A  TiON 


a  an  integration  limit 

b  »  c.,  _d,  function*  of  m. 

ninini  ni 


n£l'n£2'n£3, 
nf4’nf5'n£6’  ^ 

n h 


functions  of  m. 

n  i 


d 


water  depth 


g(m),  f(m) 
h(m) 


functions 


g 

i»  jf  h,  n 


m 


nmi 


P 

/ 

P 


r 


t 

t' 


x,  y,  s 


acceleration  of  gravity 
summation  indices 
dimensionless  wave  number 

zeros  of  the  integrand  (i.e.,  the  m.  are  the  ordered 

zeros  of  f(m)  n  1 

n 

dimensionless  pressure 
pressure 

dimensionless  radial  distance  from  the  origin 

radial  distance  from  the  origin 

distance  from  a  point  (r,  9)  to  a  aource  point  (r  ,  6  ) 

o  o 

space  variables  of  a  source  point  (cylindrical  coordinate) 

dimensionless  time 
time 

dimensionless  space  variables  (Cartesian  coordinate 
aystem) 
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\ 


x',  y\  *' 

WB 

A'  nAl'nAi' 
nAr  nAi 

B'nBl’nBi' 

nBI’„ai' 


I 

I 


C,nCl'nCi 


D*  nDl  *  nDi* 
nDl*  „Di 


■pace  variables  (Cartesian  coordinate  system) 
dimensionless  bottom  velocity 

constants 


constants 


constants 


constants 


E'nEr«Ei’ 

E/  g/  constants 

n^l’n^i 


G 

IQ 

W1* 

ln 


Jn(mr) 

R 


S 

T(t) 


time-dependent  Green's  Function 
a  function  describing  an  initial  Impulse 
integrals 

the  number  of  seros  of  the  integran 
saro  order  Hankel  transform  of  I 

o 

n**1  order  Bessel  function 

distance  between  points  (x»y*t)  and  x  , y  ,*  );  also  a 
parameter  of  the  initial  disturbance 

a  set;  also  a  surface 

a  function  describing  the  time  dependence  of  the  Initial 
deformation 
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♦ 

i 
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******** 


V 

V' 


tt 

n 


l'n 


a2*  nft3 


0 


V 


p 

T 

T* 

r' 


0 


dimensionless  velocity 

velocity 

functions  of  nn^ 

an  integration  limit 
Vm  tanh  m 

dimensionleas  wave  amplitude 
wave  amplitude 

n^1  order  Hankel  transform  of  Tl 

o 

function  describing  the  initial  deformation 

dimensionless  wave  length 

wave  length 

constant 

water  density 

a  dimensionless  time  variable  of  the  source 
a  constant 

a  time  variable  of  the  source 

velocity  potential  (subscripts  denote  partial  differentiation) 
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